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INTRODUCTION 

In developing a mathematical model for a tandem rotor helicopter 
to use in the analysis of an automatic approach and landing system for 
that helicopter, linearized equations not normally used in aircraft 
stability analysis were derived. These equations are derived herein 
along with those normally used (refs. 1,2). 


DERIVATION OF THE NON-LINEAR EQUATIONS OF MOTION 

Before developing the linearized equations of motion for analysis 
purposes, the non-linecu: equations are developed. 


Rotary Motion 

Consider the motion of an aircraft as shown in this figure. 


Motion 

Referred 



Assume that the x-z plane is one of symmetry and define the inertia 
matrix as 
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XX 
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J 


xz 
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0 



J 


xz 


0 


0 



-1- 


where 


J^y ^ / xy dm = 0 

J^y = -/zy dm = 0 

J = - f xz dm 
xz J 


due to x-z symmetry 


Note in Figure 1 that 


w = Q 


where P, Q, R is the angulcur velocity resolved along the body axes 
X, y, and z, respectively. Similarly, let the applied torque be 
denoted as 


’L" 
= M 


With the assumption that the angular momentum of any rotating machinery 
on-board, such as engines and, in helicopters, the rotor, is negligible, 
the equations in inertial space are written simply 

f = ^ (lu) . (1) 

On the assumption that the Earth is inertial, then Eq. (1) is written 
in the body axis system as 


T = I + wxlw 


where the prime denotes "as seen by an observer fixed in the body 
axes," and therefore 


P 

at “ 

-R 
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where the dot denotes tota] time derivative. Substituting the appro- 
priate definitions into Eq. (2) results in the equation set 




• 
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Jxz« 
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"xx^ 
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■> « 


.^XZ^ ^ 

I R 
zz 
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These equations define the rotary motion of the aircraft. 


(3) 


Trans latory Motion 

Assume that the mass of the aircraft is constant and denoted m. 
Let the force vector be denoted 

"X" 

_ A 
F = Y 

Z 


where X, Y, and Z, lie along the body axes. Note from Figure 1 that 
the total velocity 

U' 

V 
W 



is resolved along the body axes. The equation in inertial space may 
then be written 

F = ^ (mV^) (4) 


On the assumption that the Earth is inertial, then Eq. (4) is written 
in the body axis system as 

d'V _ _ 

F = m + wxrnV . (5) 

dt T 
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Making appropriate substitutions into Eq. (5) gives the equation set 


‘X* 


"mU" 


“mQW - mRV' 

Y 

= 

mV 

+ 

mRU - mPW 

-Z. 


_mW_ 


_mPV - mQU_ 


These equations define the translatory motion of the aircraft. 


Euler Angle Transformations 

Consider the aircraft of Figure 1. We wish to describe the 
attitude of the aircraft with respect to a set of axes fixed in the 
p^arth. To do this we define an Euler angle set denoted , Q, 

^ hese three angles are the azimuth change, elevation change, and roll 
required to arrive at the aircraft attitude from the inertial axes. 
They must be taken in the given order. If a vector is denoted in the 
inertial coordinates as C. and viewed from a coordinate system which 
has been slewed . through we get 




where 




cos 

-sin 'F 
0 


sin ^ 
cos Y 
0 


0 " 

0 ' . 
1 


If the observer is then elevated through 0, he sees 


(7) 


^2 ^y^®^^l = Ty(0)T^(’F)C^ (8) 

^here 

’cos 0 0 -sin 0" 

0 10 
_sin 0 0 cos 0, 

Finally, if the observer is rolled through $, he sees from aircraft 
body axes 


A 

Ty(0) = 
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where 




( 9 ) 





0 

cos <I> 
-sin $ 


0 

sin $ 
cos $. 


Gravity Forces 

With Eq. (9) the gravity force may be written in body axes as 


\ = T^(*)Ty(0)T^(y) 

1 

0 o 

1 

= mg 

‘-sin 

sin 

0 

$ 

cos 6 




. cos 

$ 

cos 0. 


Euler Angle Rate Equations 

The body axis rates may be written as functions of the Euler 
angle rates : 











■ 

■p" 



T 



‘o' 


‘0‘ 








• 




Q 

= T^(«) 


0 

+ T (0) 


0 

+ T^(V) 

0 






JT 




• 


_R. 



.0. 


- 

.0. 


M md 




'l 

0 

-sin 

0 

* 

_ _ 

4 

s 

0 

cos 0 

sin 

<& 

o 

o 

to 

0 

• 

0 


.0 

-sin 0 

cos 

$cos 

0. 

• 

J. 


Inertial Velocity Equations 
If the velocity in inertial coordinates is denoted 



V 

z 


(10) 


( 11 ) 
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then we may write the body axis velocities as functions of the inertial 
velocities 


u 

V 

= T^($)T„(0)T,{4') 

1 

X 

w 

X y z 

L-;j 


: 


V cos 0COS 4' + V cos Osin H' - V sin 0 
X y z 

V^(sin $sin 0cos 4' - cos $sin 4') 


V (cos $sin 0COS 4* + sin $sin 4') 
X 


+ < 


.1 


V (sin <I»sin Osin 4' + cos $cos 4') + V sin $cos 0^ (12) 

y z 

V (cos $sin Osin 4' - sin $cos 4') + V cos 0cos $ 


LINEARIZING OF EQUATIONS USING DELTA PERTURBATIONS 

The purpose of this linear model is to analyze an automatic 
approach and landing system. As such, the system is tied to the 
Approach Navigation Frame. This is an Earth-fixed coordinate system 
with the origin at the desired touchdown point, the X axis along the 
runway, and the Z axis down along the local vertical. The system is 
considered inertial. The variables to be commanded will be in the 
ANF such as Vx, Vy, and Vz and 4', 0, 4>. The- control system will 
therefore be feeding back these quantities from an inertial platform. 
The linearized equations are then desired in terms of these variables. 
The derivation of these equations follows. 


Body Trans latory Equations 

Rewrite equation set (6) with the forces being divided into 
►aerodynamic forces and the gravity forces of Eq, (10) 

^A 

U=RV-QW+— -g sin 0 
m 

V=PW-RU+-^+g sin tcos 0 (13) 

m 


W=QU-PV+— +g cos $cos 0 
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Assume the following perturbations 


U = Uq + AU P = AP 

V = AV Q = AQ 


0 = 0Q + A0 
$ = A$ 


W = Wq + AW R = AR 


= AH' 


(14) 


Note that the perturbations are simply added to the Euler angles. 
If steady flight is assumed, the steady flight equations are obtained 
from Eq. (13) : 


0 = 


^AO 

m 


- g sin 0, 


'AO 

m 


(15) 


0 = 


“AO 

m 


+ g cos 0, 


Next note that 


sin (0Q + A0) = sin 0 q + (cos 0^) A0 
cos (0Q + A0) = cos 0Q — (sin 0^) A0 


} 


(16) 


Substituting the definitions of Eq. (14) into Eqs. (13), using Eq. (16) 
and subtracting off steady flight, Eq. (15), yields the following pertur 
bation equations 


AX 

AU = -W.AQ - (g cos 0,,)A0 + — r 
0 u m 


ay, 

AV = WqAP - UqAR + (g cos 0 q)A$ + 


(17) 


AZa 

AW = UqAQ - (g sin 0 q)A 0 + -jjp 



Body Rate Equations 

Substituting the perturbation definitions of Eq. (14) into 
Eq. (3) and dropping second-order terms gives 


AP 



AR + 


AL 

I 

XX 


AQ 

AR 


AM 

^yy 

- AP 

■*■22 



( 18 ) 


Euler Angle Equations 

The relation between Euler angle perturbed rates and body per- 
turbed rates is desired. This is obtained by inverting equation set 
(11) and applying the perturbation definitions to give 


A6 = AQ 


A* = AP + (tan 0 q) AR 


AV 


AR 

cos 0Q 


(19) 


Equation Set with Aerodynamic Partial Derivatives 

The motion of the aircraft is described by Eqs. (17), (18), and 
(19). The force and moment perturbations are, in general, functions 
of the motion and can be written to first order as a linear function 
of the motion variables. Motions in the longitudinal plane are 
^assumed to separate from the lateral-directional to give the total 
^equation set which follows. 


AU = -W,.AQ - (g cos 0-)A0 + -^AU + -~AW + -^AQ + 

0 u m m m m 


AV = W.AP - UftAR + (g cos 0^)A4 + —AV + -^AP + -^AR + 

00 u mmmm 
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AW = ■ ^9 sin 0rt)A0 + — AU + —AW + -^AQ + —6 

w u m m in m 

AP = - ^AR + ^AV + y"^AP + y^AR + y-^6 

XX XX ^xx ^xx ^xx 

”u ”w “o 

AQ = ^AU + —-AW + ^AQ + ^6 

yy yy yy ^yy 

. Jv, . N„ N N. 

AR = - j^AP + j-^AV + =-£aP + y"^AR + =-^6 

22 Z2 ^Z2 ^22 ^22 

A§ = AQ 


> (20) 


A$ = AP + (tan 0 q)AR 


AV = 


AR 

cos 0, 


Equation Set with Body Rates Eliminated 

The equation set (20) may be reduced in variables by eliminating 
the body rates with the last three equations. The result is 


• \ *w ( *o\ • *6 

AU = -^AU + -^AW - W- - -?)A0 - (g cos 0.)A0 + 

m m \0m/ o m 

iO = + (w„ + i)ii + (g cos e„)4» - ^sin e„(w|) + 

+ <=“ °oK - ^))« * ^ 


z Z / Z \ . 

AW = -^AU + -~AW + (Un + -;?)a 0 - (g sin 0 ,) 
m. m \um/ ^ o 


A0 




At = 


L L 

=-^AV + ^At 

■*’XX ■‘’XX 




j-E-sin 0„ 

XX 


:cs 0 


XX 


o)« ^ I 


t, 

XX 


>( 21 ) 


•• **u ”w “a • 

A0 = j-=-AU + j^AW + J-3-A0 + -2-^ 

yy yy yy ^yy 


AY = 


®0 ■ -0 


67 (V'^ - "’xz‘* ^ V* 


-(N sin 0Q - N^cos eQ)Ay + N^s) 


in 


Reduced Equation Set in Vector -Matrix Notation 
Putting the equation set (21) in vector-matrix notation results 


Longitudinal Equations 
(ANF co-ordinate system) 
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•» ■ 



( 22 ) 
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Lateral-Directional Equations 
(ANF co-ordinate system) 
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where B = I_ cos 0- - J sin 0« and S ^ ^ the Laplace differential 

ZZ U XZ ii Qw 

operator* 


Inertial Velocity Equations 

If the equation set (12) is linearized, the result is given 


AO = AV^cos 0Q - AV^sin 0^ - WqA0 


1 


AV = AVy + WqA$ - (OqCOS 0q + WpSin 0 q)A'P • ^ (24) 


AW = AV^sin 0Q + AV^cos 0 q + UqA0 


With these equations, the body axes velocities can be eliminated from 
the equation set (21) and the result is given 


AV„ = 


- (^0 - ^0 * ® ® o )“ * 


COS 


AV„ = 


m- y m- ^ ® ®o)‘* ’ ®0 


-AV„ + 
Y. 


^■A4 + {' 
m \ 

- ®o)''» - ®o * "o’^” ®o)” 


6 

m 


11 - 



AV_ = 


A0 = 


-(tan 8„)AV^ + + ^tan 6„)4V^ - (^tan 0(, - 

* ■ (■;?"n ■ "ST^o * 9 ®n)*® * Tif*) 

COS 0Q\in VmO mO^ 0/ ***/ 

/M M V /^M \ 

(^os 0Q + j^in Go)av^ - (^in Oq - j^os G^) AV 

V yy yy / \ yy yy / 

M /M M \ M- 

+ ^3-40 - (j^o - r^o)''® * I-' 
vv V yy yy / yy 


>( 25 ) 


.. L L L 

A* = + =-2-A$ + j-^oA$ 

■"^XX ^ ■‘■xx ■“'XX 


■ ®o) 


A¥ 


XX 
L 


- 0 - ^os 0 ] A¥ - cos 0j 

\ XX XX / XX ' 

T. 

+ WjjSin 0 q), 


lAf + 


Af 


t = Kv^ - •’xz''* * V* * ®o 

- N^cos 0 q)AY - N^(UqCOs 0q + WQsin 0 q)AT + 


where again B = I „cos 0. - J sin 0^ . 
^ zz 0 xz 0 


Inertial Velocity Equations in Vector-Matrix Notation 

The equation set (25) may be put into vector-matrix notation 
with the result 
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Longitudinal Equations 
tANF coordinate system) 



Lateral Equations 
(ANP coordinate system) 



where - U^cos 0 q + W^sin 0 q, B - ®0 " ®0 ® 

differential operator. 


(26) 


(27) 



LINEARIZATION OF EQUATIONS USING EULER PERTURBATIONS 

In the previous section linearization was done to develop a 
math model for approach and landing system analysis. As such, it was 
desired referenced to the ANF. In this section the equations are 
developed for use in the analysis of a stability augmentation system. 
Such a system will generally feed back only body-referenced rates 
P, Q, R and U, V, W. The linearizations then are desired in terms of 
these variables. 

In the previous section linearization was done by representing 
the perturbed angles as the initial angles plus delta perturbations of 
the Euler angles so that 

0 = 0Q + A0 

$ s= 

¥ = 'I'q + A'i' 

In this section the perturbed angles will be arrived at by taking the 
initial angles and adding to them another set of Euler angles refer- 
enced to the initial body axes. These perturbation Euler angles will 
be denoted \|^, 6, ()>. The translational velocity and angular velocity 
perturbations will be as in Eq. (14). 



Body Trans latory Equations 

The perturbed equations will differ only in the gravity force 
which is described with , 0, $. This perturbed force is written 




‘-sin ©q‘ 


‘-sin 0 q‘ 

AF 

grav 

Tx(4>)Ty(6)Tz(i|^) 

0 

- 

0 



cos Qf. 


cos 0rt 



L OJ 


L 0 J 


" Carrying out the matrix multiplication we get 


AF 


grav 


-0COS 0Q 

i|;sin 0Q + <|>cos 0 q 
- 0sin 0Q 


mg 


(29) 


(30) 


- 14 - 



The perturbation equations then become 


AU 

AV 

AW 


-WqAQ - (g cos 0 q)6 + 

WqAP - UqAR + (9 cos 0q)4' + (9 sin 0q)i|' + -jjj— 

AZ 

UqAQ - (g sin 0 q)6 + -~ 


( 31 ) 


Body Rate Equations 

Since angles are not involved, these equations are just those of 
set (18) . 


Euler Angle Equations 

With the defined transformations, the body rate perturbations 
can be written 


AP" 

AQ 

= T^(*) 


1 

0 -©-. 
1 

+ T (6) 

y 

1 

0 •CD 

+ 

'o' 

0 

m 

AR. 



LoJ 


LoJ 


i- 

. 


With small angle assumptions on the perturbations and by dropping 
second-order terms the result is 


AP" 


Y 

AQ 

= 

• 

0 

AR. 


i. 


(33) 


This means that for this linearization, unlike that of the prior 
section (see Eq. (19)), the angular perturbations are simply the 
integrals of the body rates. Although, for convenience in the remain 
der of this section, the angles ijj, 0, $ will be used, their simple 
relations to the body rates should be remembered. 


Equation Set with Aerodynamic Partial Derivatives 

Using assumptions similar to those of set (20) and Eqs. (33) to 
eliminate body rates yields the following equations: 
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X X X X 

AU = -W^e - (g cos 0-)e + ~AU + —AW + -^AQ + —6 
u u m in in m 


AV = Wq^ - Uq^ + (g sin 0 q)iJ) + (g cos 0 q)(J) + -j^AV + 


+ + “6 
in m 


. z z z z. 

AW = U_e - (g sin 0»)0 + — AU + -^AW + Aq + —5 
u u in m m m 




^ ^ + ^Av 


XX 


XX XX 


XX 


XX 


(34) 


M . M M M. 

0 = =-^ + j-Si-AU + ^AW + j-^6 

yy yy yy yy 


•^xz- ®^p* ”r- ** 

^ + _£h(. + + ^AV + - 


zz 


zz 


zz 


zz 


_6_^ 

zz 


Reduced Equation Set in Vector-Matrix Notation 

Putting the equation set (34) in vector-matrix notation with 
Laplace notation results in a set of longit’idinal equations identical 
to set (22) for the delta perturbations. 


Longitudinal Equations 
(body axis system) 
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AW 


0 

\ 0 m/ ^ 0 




m 

/ M \ 




Mjt 

(s-^]s 


0 



\ I / 




I 

1 

>1 

>1 




yy 
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The lateral equations, however, are different. 


Lateral-Directional Equations 
(body axis system) 





sin 0, 






XX 


6 (36) 


where 


^ dt ' 


the Laplace differential operator. 


CONCLUDING REMARKS 

Two sets of linearized equations of motion for an aircraft 
initially in steady flight have been derived. The first set is written 
in terms of Euler angles and approach navigation frame velocities for 
the purpose of analyzing an automatic approach and landing system. 

The second set is written in terms of. body rates and body-referenced 
velocities for the purpose of analyzing simple stability augmentation 
systems. The two equation sets are clearly different.. The difference 
depends on the aircraft and its steady flight condition. Any assump- 
tion that the second more common set may be used in place of the first 
is not, therefore, automatically correct for all aircraft and all 
steady flight conditions but should be verified in each instance. 


REFERENCES 


1. Tandem Helicopter Trim and Stability Analysis, Boeing-Vertol 
Division Report Aero. Inv. III-264, March 29, 1965. 

2. Dynamics of the Airframe, BuAer Report AE-61-4II, September 1952. 


- 17 - 



